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THE CRACK-CONTACT AND THE FREE-END PROBLEM 
FOR A STRIP UNDER RESIDUAL STRESS (*) 

by 

M. BaJcioglu and F. Erdogan 
Lehigh University, Bethlehem, PA, 18015 


ABSTRACT* The plane problem for an infinite strip with two edge cracKs under 
a given state of residual stress is considered. The residual stress is com- 
pressive near and at the siarfaces and tensile in the interior of the strip. 

If the crack is deep enough to penetrate into the tensile zone, then the prob- 
lem is one of crack-contact problem in which the depth of the contact area 
is an unknown which depends on the crack depth and the residual stress profile. 
The problem has applications to the static fatigue of glass plates and is 
solved for three typical residual stress profiles. In the limiting case of 
the crack crossing the entire plate thickness, the problem becomes a stress- 
free end problem for a semi-infinite strip under a given residual stress state 
away from the end. This is a typical stress diffusion problem in which decay 
behavior of the residual stress near and the nature of the normal displacement 
at the end of the semi-infinite strip are of special interest. For two typ- 
ical residual stress states the solution is obtained and some numerical re- 
sults are given. 


1. Introduction 

Introducing residual stresses into structural components which are 
compressive near and at the surfaces to improve their impact and fatigue 
resistance has been a design practice for many years. Some of the processes 
used for this purpose are tempering, cladding, ion exchange, and shot peen- 
ing. In calciilating the stress state in the part, these residual stresses 
must be superimxjosed on the stress state resulting from the applied loads. 

In some cases residual stresses may be the only stress state in the body. 

For example, in considering the problem of subcritical crack growth due to 
static fatigue in glass plates and other ceramics which normally do not 
carry any external loads, the crack driving force is mainly provided by the 
residual stresses. Such a subcritical crack propagation may take place if 
the siirface crack accidentally introduced into the plate is deep enough for 
the initial crack front to be in the tensile zone and if the corresponding 
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stress Intensity factor is greater than the threshold level which is the 
ainimum stress Intensity level required for crack growth under static 
loading (for a review of the subject see [l]) . The correlation between 
the subcritical crack growth velocity and the stress intensity factor in 
ceramics, polymers, and certain metallic alloys under adverse environmental 
conditions appears to be quite well established [l,2]. Therefore, in prin- 
ciple it is possible to make a reliable prediction for the failure or crack 
arrest time in the material under residual stresses provided the proper 
fracture mechanics analysis is available. 

In this paper a relatively simple problem of a plate xmder a known 
state of residual stress is considered. The plate is assumed to have two 
symmetric edge cracks. However, because of the compressive stresses, near 
and at the boundaries, the crack surfaces will be closed along a certain 
unknown distance from the boundary. Therefore, the problem is one of crack- 
contact problem rather than a conventional crack problem. After the cracks 
go through the entire plate thickness, the problem reduces to a stress-free 
end problem for a semi-infinite rect. ngular strip under residual stresses. 
This is a typical stress diffusion problem and is treated as the limiting 
case of the edge crack problem. To study the behavior of the displacement 
at the stress-free end and of the diffusion of residual stresses is one of 
the primary aims of this paper. The plane problem of an infinite strip 
with various crack geometries has been considered by many investigators who 
used a variety of techniques to solve the problem [e.g., 3-8], 

2. On the Formulation and Solution of the Problem 

Consider the elastostatic plane strain or generalized plane stress 
problem for an infinite strip under a symmetric state of residual stress 
satisfying 

dyy(x,y) = a^(x) = a^[-n) , I aj^(x)dx = 0 , (1) 

-h 

where d_(x) is a known function (Figure 1). Let the strip contain two 
R 

symmetrically located edge cracks along y=0, a<|x|<h. Since the stress 
component perr^endicular to the crack is compressive near and at the 

surfaces, the crack faces will be closed along b<|xl<h, where the constant 
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b is not ycnown. The contact of the crack surfaces at |x| =b is "smooth" 
and, consequently, the "ends" of the slits at |x[ =b have a cusp shape 
rather than the standard parabolic form. Note that the stress intensity 
factors at the ends of the slits shovm in Figure lb are defined by 

= IJj [2(a-x)]‘’clyy(x,0) - ^ jji» liCk-al]'' ^ ' 

■S - iiS [2(x-b)]’’o^y(x,0) - - ^ Um [2(b-x)]'’ Ifvlx.o) (2X.W 

where v is the displacement component in y direction, and p and < are the 
elastic constants, K= 3-4v for plane strain and K= (3-'V)/(l+v) for general- 
ized plane stress, V being the Poisson’s ratio. Thus, from (2b) it may be 
seen that the condition of "smooth closure" at the end of the slit x=b may 
be expressed as 



(3) 


Equation (3) provides the additional condition to determine the unknown con- 
stant b. 


The solution of the edge crack problem may be expressed as the sum of 

two solutions! A) the homogeneous solution in the strip without any cracks 

which is essentially given by (1) , and B) the perturbation solution for the 

B 

strip with the edge cracks in which the crack surface traction 0yy(x,0) = 

- 0 (.x) is the only external load. It is clear that problem B contains all 
the important information for the stress-free end as well as the crack prob- 
lem. The formulation of the problem is identical to that given in [s] where 
it was shown that the problem described in Figure lb may be formulated in 
terms of the following integral equation which, considering the symmetry, is 
written for one quarter (y>0, 0<x<h) of the medium only; 


4y 

TT(1+IC) 


,b 

„‘-t-3 


+ k(x,t)] G(t)dt = cJ® (x,0) =- 


t+x 


yy 


C^tx) 


a<x<b 


(4) 


where 


G(x) = v(x, +0) , a<x<b , 


k(x,t) = 


[K(x,t,s)e^® - K(x,-t,s)e ^^]ds , 


-ts- 


(5) 

( 6 ) 
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K{x,t,s) - {[hs-X2hs-3) (h-t) s-(2+ts) e ^^®3cosh(xs) 


“ [l-2s(h-t)+e ^^®jxs sinh(xs) }[2hs+sinh(2hs)] , 


( 7 ) 


In the general case the index of the singular Integral equation is +1 and 
hence its solution is determinate within an arbitrary constant [9], On the 
other hand, in deriving (4) it was assumed that for 0<[x|<a and b<lx|<h G{x) 
rather than v(x,0) vanishes. Thus from (5) it follows that 


G(x) dx = 0 


( 8 ) 


which is used to determine the constant arising from the solution of (4) . 

The solution of the contact-crack problem is then obtained by evaluating 
G(x) and b from (4) and (3) subject to the verification that the resultant 
stress 0^y(x,0) obtained from the superposition of the solutions of problems 
A and B for b<lx|<h be compressive. 

It should be noted that the left hand side of (4) gives the normal 
stress a within as well as outside the crack on y=0. Therefore, after 


yy 


obtaining G(t) and b, (4) may be used to evaluate o^y{x,0) on 0<|x|<a and 
b<|x|<h. In the residvial stress problem for a>0 b is always less than h, 
and b^h for a-+0. This limiting case of a=0, b~h corresponds to the "end 
problem" for which (4) is still valid. However, in this case the kernel 
k{x,t) contains additional singular parts which must be properly separated. 
The technique for doing so was described in [sj and will be omitt ' • ' this 
paper. The expression which will be needed to study the stress ditfc':>ion 
phenomenon is that of a^^(x,y), 0<x<h, y>0. From the general formulation 
of the problem it may be shown that [5,7] 

,b 


V ^ TfO 

4y yy 


(x,y) = 2 G{t)dt (l+ys)e ^®cos xs sin ts ds 

J A ^ r\ 


G(t)dt 


[K(X,t,! 

JO 


s) e^^ - K(x,-t,s) e '^^Icos ys ds (9) 


-ts- 


where K(x, t,s) is given by (7) . For y=0 and a<x<b (9) reduces to (4) . 
Evaluating the inner integrals, the kernel in the first term on the right- 
hand side of (9) may be expressed in closed form and the first term becomes 
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t-x 

(t-x) 


+ 


t+x 

(t+X) 


4 



For the limiting case b=h the kernel in the second term of (9) becomes 
unbounded as t 2 uid x approach h. For correct numerical land singularity) 
analysis this unbounded part needs to be separated. After doing so, for 
b<°h (9) may be expressed as follows: 


14< 

4b 



(x,y) 


’^r t-x t4-X 2y^(t-x) 

^(t~x)"^4-y^' Tt+xP^+y^ [ (t-x) '^+yj^ 

• A 


4 . 2y^ (t4X) 2(2h-x-t) (4h-x-3t) [ (2h-x-t) ^-y^! 

nt+xp+yT” “ (2h-x-t)'=+y^ L (2h-x-t) ^+y^ 


4 4 (h~t) (x-h) 


(2h-x-t) ^-ly^ (2h-x-t) 

L{2h-x-t)ny^P 


G { tl dt 


rh f«> ts 

G{t)dt {[K(x,t,s)-K (x,t,s)]e 
J a jo ®° 


- K(x,-t,s)e }cos ys ds , (0<x<h, y>0) (10) 

where K (x,t,s) is the asymptotic behavior of K(x,t,s) and may be obtained 
from (7) as 


K^(x,t,s) 


[-2+s{4h-x-3t) 4 2s^ (h-t) (x-h) ]e ® 


( 11 ) 


Another quantity of interest is the end displacement in the semi- 
infinite strip which referring to (5) may be evaluated from 


v(x,0) - v(0,0) 


G(t)dt . 
. o 


( 12 ) 


For 0<a<b<h the numerical solution of (4) is straightforward, although an 
interpolation scheme is needed to determine the unknown constant b from the 
condition (3). For a-^O, tr+h and the crack problem becomes a stress-free 
end problem. In this case G(t) may be evaluated by using a numerical tech- 
nique similar to that described in £4]. 
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4. Results for the Crack Problem 

The solution of the problem is carried out for three different symmetric 

(•) 

residual stress distributions given by 


ohx) = 0^(l-3xVh^) , 

K O 

(13) 

<J^(x) = a^d-SxVh**) , 

(14) 

o^(x) « a (i-7xVh®) , 

K O 

(15) 


where a is the magnitude of the tensile stress in the midplane x=0 (Fig- 
o 

vure 1) . First the problem is solved routinely only for those values of a 

and b for which the stress intensity factors and k^^ obtained from (2) 

are both positive. The results corresponding to the residual stresses (13) , 

(14), and (15) are given in Figures 2, 3, and 4, respectively. The stress 

intensity factors shown in these figures are normalized with respect to 

where l-(b-a)/2. The figures also show the tensile portion of the 

numerical stress 0^/0 , (i=l,2,3). Note that 
R o 




as a-^b , (i=l,2,3) 


This is due to the fact that for small values of S. = (b-a)/2 the problem is 
equivalent to that of an infinite plane containing a crack of length 2£ and 
siibjected to crack surface traction -0^ for which k =)t =o (i=l,2,3). 

K a D H 


For a given value of a, the value of b for which kj^=0 (i.e., the loca- 
tion of crack closure and beginning of contact area) is shown in Figure 5. 
From the figure it may be observed that the point a=b on these ctirves corre- 
sponds to the point at which the residual stress is zero, 

(i=l,2,3), (Figure la). It may also be observed that the value of b for 
which kj^=0 is always greater than x^. For a given value of the crack depth 
a the stress intensity factor corresponding to kgO is shown in Figure 6. 
{*) 

The parabolic distribution (13) seems to be typical for residual stresses 
in tempered glass and the 6th degree polynomial (15) more representative of 
the internal stress induced by ion exchange in glass plates. 
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!lote that as expected k -+«' as a-*-0 and k -►0 as a-oc . 

a a o 

Figure 7 shows a sample result giving the stress distribution 0yy(x,0) 
on y=0 plane for the crack-»contact probleir As pointed out before^ the 
total stress is 


where 

0^^(x,0) = o^(x) , (18) 

The particular residual stress lised in this example is given by (14) and is 

also shown in the figure. Since O’ (b(0)=0, the solution of the perturbation 

B 

problem must give a (b,0)--a^(b) . This was found to be the case in the nu- 
merical solution within four significant digits. 

4. Results for the End Problem 

In the end problem the main interest is in the diffusion of residual 

stress o (x,y) in y direction# going from o =0 for y=0 to 0 =0 (x) for 

yy yy yy R 

y*". This is shown in Figures 8 and 9 for the residual stresses given by 

(13) and (15) , respectively. These results as well as those given in Fig- 
ures 10 and 11 must be considered in conjunction with (17) and (18) . Fig- 
ures 8 and 9 show 0^{x,y) for y/h=0, 0.4, 0.8, 1.2. Perhaps a better 
description of the stress diffusion may be observed from Figures 10 eind 11 
where 0^ at x=0 and x=h is given as a function of y. Examining the figures 
one could make the following general remark; the stress diffusion rate on 
the surface is greater than that in the mid-plane, and a plate thickness 

away from the end, i.e., at y=2h the stress at the surface drops to approx- 

{*) 

imately 1% of its maximum value (which is at y=0) . 

To give an idea about the deformed shape of the stress-free end. Fig- 
ure 12 shows some sample results. Here the relative displacement v(x,0) -v(0,0) 

(*) 

Needless to say this depends on the undisturbed residual stress profile 
Oj^(x) . For example, for Oj^ given by (13) this figure is approximately 1.1% 
and for given by (15) it is 0.83%. 


- 7 - 



is gxven for residual stresses (13) and (Id) > The normal ixatxon factor 
which appears in the figure is 


(l+K)ho 

V T ^ 

o 4y 


( 19 ) 


The examination of the results shown in Figures 8-12 appears to indicate 
that in relative terms the difference between two end displacements is 
much greater than that between two corresponding residual stresses. 
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Stress intensity factors at the crack tips a and b due to th'.- 
residual stress field Og(x) * c (l-Bx'/h" ) . 


Figure 3. 
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Figure 7. The normal stress on the y*0 plane for the crack-contact probl 
crack: a<x<b, contact: b<x<h. 
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Figure 11. Variation of Oy^ on x=0 and X“h in a semi-infinite strip under 
the residual stress o^. 
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Figure 12, Relative normal displacement of the end surface in a soni-inf mite 

strip under the residual stresses and (v *(lfic)ho /(4p)). 
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